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Abstract 

We perform calculations for the production of (j) mesons in nuclei at energies just 
above threshold and study the A dependence of the cross section. We use results 
for the (j) selfenergy in the medium obtained within a chiral unitary approach. We 
find a strong A dependence which is tied to the distortion of the incident proton and 
to the absorption of the (j) in its way out of the nucleus. The effect of this latter 
process reduces the cross section in about a factor two in heavy nuclei proving that 
the A dependence of the cross section bears valuable information on the (p width in 
the nuclear medium. Calculations are done for energies which are accessible in an 
experimental facility like COSY. 

1 Introduction 



The study of the properties of vector mesons in a nuclear medium is one of the subjects in 
hadron physics which receives continuous attention (see for instance Ref. [1]). Although 
originally the p meson properties were mostly investigated, with the time the renormaliza- 
tion of the properties has been taking over as it is linked to the way kaons and antikaons 
are modified in the nuclear medium, which has also been the subject of intense study. One 
of the motivations for this latter work is the deviation of the K~ selfenergy from the low 
density theorem, which is needed to reproduce kaonic atoms data [2-5], together with the 
possibility of formation of kaon condensates in neutron stars [6] . 

Another reason why the (p properties have got a renewed interest is the fact that the 
medium renormalization in the case of the is relatively more drastic than that of the p. 
Indeed, predictions of an increase of the width by a factor five or six [7, 8] to ten [9] , 
at normal nuclear matter density, have been made using different chiral approaches. This 
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large width should in principle be detectable experimentally and, indeed, different reactions 
have been studied or suggested like (f) production in A — A collisions [10, 11], the reaction 
7r~p — > (pn in nuclei [9] or different methods based on inclusive (j) photoproduction in 
nuclei [12-14]. 

The width of the has not been the only matter of concern and the possible shift of 
the mass has also captured attention. In that sense, the mass change has been studied 
in several approaches like using effective Lagrangians [9,15-18], QCD sum rules [19,20] 
or the Nambu-Jona-Lasinio model [21]. In ultra-relativistic heavy ion collisions the (j) 
width modification in matter has also been studied in a dropping meson mass scenario 
[17,22-26], as a result of coUisional broadening through 0-baryon [27] or 0-meson [28] 
scattering processes. All these works point at a sizeable renormalization of the width 
and a small mass shift^. For this reason we will only be concerned about the width in 
the medium. 

The aim of the present work is to propose a method to determine the width in 
the nuclear medium. The traditional method in the works quoted above (except [14]) 
is to look for a broadening of the width reconstructed from the invariant mass of its 
decay products. Here, instead, we use a different philosophy and we investigate the A 
dependence of (j) production in pA collisions, in a similar way as it was done in [14] with 
the (j) photoproduction in nuclei, which is the subject of experimental investigation at 
Spring8/0saka [29]. The advantage of performing the reaction slightly above threshold is 
that one can rule out the contribution from coherent production which might obscure the 
interpretation of the experimental results in [29]. The present reaction, with its particular 
kinematics, is amenable of experimental performance at facilities like COSY. 

2 Nuclear effects in the (j) production 

2.1 General formalism 

In order to implement the relevant nuclear effects in the production cross section we 
will use a model based on many body techniques, successfully used in the past in many 
works [30,31] to study the interaction of different particles with nuclei. The model assumes 
a local Fermi sea at each point in the nucleus and provides a very simple and accurate way 
to account for the Fermi motion of the initial nucleon and the Pauli blocking of the final 
ones. On the other hand, we have to take into account the distortion of the incoming 
nucleon and the final meson in the their way through the nucleus, which are evaluated in 
the present work using an eikonal approximation. In fact, the A dependence of the effect 
due to the final absorption will provide the way to test the modification of the meson 
width in nuclear matter, as we will explain in much more detail in the following, and which 
is the main aim of the present work. 

^Exceptionally in Ref. [19] a mass shift of a few hundred MeV is reported under hot and dense matter 
conditions. 
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Within the local Fermi sea approach the nuclear cross section can be evaluated, as a 
first approximation, as: 



/f d^k 
d^r J j^-^Qikp - \k\)amiPLab,k,r} (1) 

since 4 J (|^0(^f— |^|) = pi^"), where k is the momentum of the nucleon in the Fermi sea, 
kp the Fermi momentum at the local point, G the step function, piab the momentum of 
the incident proton and 0"^ the elementary pp ppcf) and pn pncf) average cross section 
in the nuclear medium, which will be defined later in Eq. 

At this point we can add in Eq. the following eikonal factor to account for the 
distortion of the incident proton in its way through the nucleus till the reaction point: 



exp 



<ypN {PLab ) P ( Vb^ + z'^) dz 



(2) 



where z and h are the position in the beam axis and the impact parameter, respectively, 
of the production point r of Eq. (|T)). In Eq. Q is the total pp and pn averaged 
experimental cross section, taken from [32], for a given incident proton momentum. Eq. Q 
represents the probability for the proton to reach the reaction point without having a 
collision with the nucleons, since (Jp^p is the probability of proton collisions per unit 
length. The use of this eikonal factor will select the one-step processes and neglect the 
possibility that the is produced in a second collision of the nucleon, or a possible excited 
A. We shall also take these two-step processes into account later on, but we advance that 
the A dependence of the production cross section is already given quite accurately by 
the one-step process. 

The final absorption in its way out of the nucleus can also be accounted for by means 
of a similar eikonal factor and for the evaluation of the probability of loss of flux per unit 
length we can proceed as follows: let n(p(^,p(r)) be the selfenergy in a nuclear medium 
as a function of its momentum, p^, and the nuclear density, p{r). We have 

^ = ^opt = ^eVopt + iXmV^opt , (3) 

and hence 

2 luJif) uj^ at 

where a;^ is the energy and V is the decay probability, including nuclear quasielastic 
and absorption channels. However, in what follows we shall only include in TmU the 
absorption channels of the 0, since in quasielastic collisions the nucleus will be excited 
but the will still be there to be observed. Hence, we have for the probability of loss of 
flux per unit length 

dV _ dV _ dV _ JmH 

~dI~Vdt~^m~ 'pT ' ^' 
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and the corresponding survival probability is given by 



exp 



f°° (-1) 
Jo 



(6) 



where r ' = r + Z-j^ with r the (p production point inside the nucleus. The study of the A 
dependence of the total nuclear cross section due to the absorption effect, Eq. ©, is the 
main aim of this work, since it reflects the modification of the (p meson width in nuclear 
matter. 

The study of the A dependence was also done for rj photoproduction in [33], along 
the same lines as discussed above, in order to determine the inelastic t]N cross section in 
a nuclear medium. That paper also shows that an alternative BUU approach [34] gives 
similar results to the Glauber, or eikonal, approach. 

Another relevant nuclear effect of sizeable consequences is the binding energy. When 
studying nuclear reactions with other beams, for instance photons, one usually neglects the 
binding energy, Vn, because it affects both the initial and outgoing nucleons and then it 
cancels in the 5 function of energy conservation. Therefore, one usually considers only the 
kinetic energy of the nucleons. In the present reaction, with a proton as incident beam, the 
same cancellation happens, because we have two initial and two final nucleons. However, 
in order to neglect Vn of the nucleons we must consider the kinetic energy of the nucleon 
beam at the reaction point which is bigger than the asymptotic value due to the attractive 
potential felt by the proton inside the nucleus. This increase in the kinetic energy can be 
evaluated considering that Etot = Ekin{asymptotic) = Ekinilocal) + Vsilocal), where Vg is 
the potential due to the local Fermi sea, which in the Fermi model is Vs{r) = —eF{r) = 

—^ij^, with kp = (|7r^p(r))-'^/^ and M the nucleon mass. With these considerations we can 
define the local initial proton momentum p^^^ such that 

kUr) 



+ P'L, = yjM^ + pl^, + (7) 

This p'l^ah is incident proton momentum used as input to evaluate the elementary cross 
section. 

The elementary cross section in the nuclear medium for pip'i^ab) + N{k) —>■ N{pi) + 
p{p2) + (p{p^) reaction is 



X E E - kFir))eip2 - A;^(r))e(l - A^) (8) 

— * 

where P = p^^^ + k — p^, and A, providing the cosinus of the angle between P and pi, A = 
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y + 5*^. In Eq. (jHJ the azimuthal angle of pi with respect to P has aheady been inte- 
grated, assuming that |Tp does not depend on this angle. This is, however, supported by 
the experiment [35] where the angular dependence of pp ppip is almost flat. Hence we 
assume in what follows |Tp to be angular independent. 

Gathering all these results the final expression for the production cross section in 
nuclei reads: 



(^Aip 



M3 



Lab 



(27r)7 \pL. 



ab 



d^k / dp 



dVtri 



d'b 



dz exp 



-WITI 



{pLab)p{VWT^)dz' 



exp 



\P\E{p,)uj{p^) 
eikp - \k\)e{p, - kF{r))eip2 - A;p(r))e(l - 



dl—^TmIl{p^, p{r')) 



(9) 



As argued above, the T matrix can be considered a constant for a given energy, which 
means that we can divide Eq. Q by the free reaction cross section to get rid of |Tp. The 
free reaction cross section can be easily evaluated as 



M 1 / 

afreeiPLab) = T^TTTTS / dE{p,)dE{p,)Q{l - B^)Q{^s - E{p,) - E{p,))J2 E 1^1' 

\PLab\ l^^J J 

(10) 

where 5 = (^i - E{p^) - E{p2)f - - p^^ - P2^ / {2\p^\\p2\) . 

The observable we will evaluate in the results section is R = a a/ (Aajree) from Eqs. © 
and (fTUI). 

For the evaluation of ImU in nuclear matter we use the results of the model of Ref. [8] 
and its extension to finite 0- meson momentum done in [14]. This model is based on the 
modification of the KK decay channel in the medium by means of a careful treatment of 
the in medium antikaon selfenergies [7]. It uses a self consistent coupled channel unitary 
calculation, based on effective chiral Lagrangians, and taking into account Pauli block- 
ing, pion selfenergies and mean-field potentials of the baryons (for the S-wave part) and 
hyperon-hole excitations (for the P-wave part). 

In Fig. ^ we show the imaginary part of the selfenergy as a function of the mo- 
mentum for different nuclear densities. The contribution to 11 coming from the free 
decay into KK, non density dependent, has been subtracted from the full 11 since the KK 
coming from the free decay would be detected and counted as a event, hence it does 
not contribute to the loss of flux required in the argument of the exponential in Eq. ©. 
As shown in [8], this corresponds to a medium width at rest at p = po of the order of 
24 MeV. 

At this point it is worth mentioning that we have used a local (momentum independent) 
nucleon selfenergy in Eq. (|7j). In practice it also has a momentum dependence which is 
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Figure 1: Imaginary part of the in medium selfenergy, without the inclusion of the free 
part, as a function of the momentum of the and the nuclear density. (From Ref. [14]). 
For a at rest and p = po this corresponds to a in the medium width correction of 
24 MeV. 



often parametrized in terms of an effective mass [36]. We propose the following expression 
for the nucleon energy in the medium 



E{p) = ^JM*'+P + M-M*-^, (11) 
which at low momenta gives the nonrelativistic expression 

E(p)=M + ^ ^ (12) 

2M* 2M* ^ ' 

used in many body approaches [36] with the Thomas Fermi energy as the local potential. 
Our expression (lllj). at the same time, provides a good relativistic energy of the nucleon 
in the low density limit. The effective mass M* is evaluated in accurate many body 
calculations [37,38], which provide a p dependence which can be approximated by 

M* = M(l-0.2p/po) • (13) 

At high energies, however, the use of Eq. (jllj) with the effective mass of Eq. (jl3p would 
overestimate the effect of nonlocalities, so our calculation using this approximation should 
be considered as an upper bound. 
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We have evaluated the production cross section modifying Eq. (|7j) to 

\Jm*' + p 'L, + M - M* - ^ = + p , (14) 

and the expression for A (after Eq. (jH))) by 

A ^ [m*' + - [E{p',J + £;(fc) - i?(pi) - uj{p,)f] , (15) 

2|P||pi| I- 

with E{p) = \/ M*^ + We shall comment in the results section that there are some 
small corrections coming from this source, but they do not modify the A dependence of 
the results. 



2.2 Two-step process with nucleon intermediate state 

Now we want to improve the former result taking into account production from two-step 
processes. Let us imagine we have a pN collision of the initial proton going to any other 
channel than production. In such cases the fast incoming proton will usually survive 
although with a reduced energy, by means of which it still can contribute to production. 
We estimate the contribution from this mechanism. The first step is to estimate the energy 
loss. The total pN cross section for Tp = 2.5 to 3 GeV is around 40 mb and consist of 
30 % elastic and 70 % inelastic, going mostly to several pion production. Estimating that 
the elastic collisions go in average around angles of = 30° or bigger in the CM frame 
since B G [0°, 30°] accounts for only 7 % of the phase space, the p energy loss in the 
lab frame is around 180 MeV or more. On the other hand, with an average three pions 
produced per collision, the amount of 500 MeV energy loss is a realistic figure. With the 
percentages given above for the elastic and inelastic collisions, one comes out with around 
I\E ^ 400 MeV loss per collision or more. 

The probability that a is produced in a second step, after a prior pN collision, is 
easily calculated assuming that the proton goes still forward after the pN collision, which 
is essentially the case in the lab frame. 

Since a^^pih^ z")dz" is the probability of a pN collision in dz" then the formula to 
evaluate the cross section for this process is given by Eq. Q substituting 



1 



\pLab\ 
1 



\Pl, 



ab 



dz exp 



'OO 

oo 



O-pTV (P Lab) pi V&^+i^) dz' 



Fip'i 



Lab) 



dz" I dz exp 



exp 



(^pN{PLab)p{ Vp+i^) dz' 



- I ^pN{PLab)p{Vb^ + Z''^)dz' 



(^pNiPLab)piVb^ + Z""^) 



Fipl 



Lab) 



(16) 



where we have introduced a new integration over z", the point of primary pN collision 
and the two exponential factors account for the probability that the proton reaches the 



7 



point z" without any collision, times the probability that the struck proton reaches the 
production point without any other collision. In this sense we guarantee that there is one 
and only one primary "pN collision. In Eq. (|T6|) the function F corresponds to the rest of 
the integral and factors in Eq. as a function of the proton momentum in the nucleus. 
Hence, v'hahi which will be the proton momentum after the primary collision is given by 



2.3 Two-step processes with A intermediate states 

In production processes below or close to threshold, multistep processes are usually impor- 
tant as we have already mentioned. Also it is known that A production in intermediate 
states is sometimes an efficient way to produce the final particles. Two methods among 
many, have particularly been used with success to deal with these multistep processes, and 
not necessarily around threshold, but particularly above it where some of the simplifications 
done are more accurate. One of the methods used is the transport equations [34,39,40] 
(BUU) in which there are collisions of the incoming particles which produce certain final 
states. In particular, A states are formed and allowed to propagate till they collide with 
other nucleons. These A's are assumed to be elementary particles, although the finite 
lifetime is taken into account. Another method which has been used in a large variety of 
physical problems is a computer simulation of the reactions [30, 41-43] , with similarities to 
the cascade models. This method, originally developed to deal with pionic reactions [30], 
has also been used to deal with pion [41] and nucleon emission [42] in photonuclear re- 
actions and in nucleon emission in electron nucleus scattering [43]. There is an essential 
difference between the transport method and the computer simulation of [30,41-43]. In the 
latter one a quantum mechanical many body treatment of the process is done in which the 
A's are never elementary particles, but they appear as propagators in Feynman many body 
diagrams, which would qualify as two or three body processes in the transport method. 
The difference is more than technical since this procedure allows one to sum over the spin 
of the A's in amplitudes (in the propagator) while in the transport method the sums over 
spins are done on the cross sections. This leads sometimes to differences in angular dis- 
tributions, the information of which might be missed in the transport method. However, 
sometimes, corrections are done in the standard BUU equations to account for these miss- 
ing angular correlations [44]. The Monte Carlo simulation of [30,41-43] evaluates quantum 
mechanically these multistep processes for the diagrams which are called irreducible. In 
this context reducible diagrams are those that, by cutting a propagator line of a stable par- 
ticle in the diagram of an amplitude, two valid diagrams result, which can be interpreted 
as a multistep process with stable particles in the intermediate states. These processes can 
be reduced to a sequence of the more elementary ones. Their probabilities are calculated 
with the rules for irreducible diagrams, and the multistep ramification is done using a 
Monte Carlo simulation procedure, allowing the individual steps to occur according to the 
calculated probability. 




(17) 
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Let us address now with the many body techniques, using irreducible Feynman dia- 
grams, the analogous of the two-step process A^A^ NA followed by AA^ NNip in the 
transport equation. 

This two body process would benefit with respect to the one considered in the former 
subsection from the fact that the A couples more strongly to pions and vectors than the 
nucleon. For instance, one can consider a mechanism from the model of [45,46] like the 
one on Fig. |2l which, with respect to the same one with a nucleon instead of a A, would 
benefit from the factor fnNA/ I-kMN = fpNA/fpNN = 2.13 in the amplitude, hence a factor 
4.5 in the cross section. 
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Figure 2: Example of one of the possible diagrams for the AA^ — NNcf) reaction [45,46]. 

We shall distinguish between A excitation on the target and A excitation on the pro- 
jectile. These two mechanisms appear in the elementary reaction pp —>■ AN [47-49], but 
they proved to be relevant in the (^He, t) reaction in nuclei [50] and in the {a, a') reaction. 
In this latter case, only A excitation in the projectile is allowed [51], which was clearly 
seen experimentally in [52]. We describe below the two mechanisms starting from the A 
excitation on the target. 

We now take the diagram shown in Fig. IHl b), which can be interpreted as having 
a A^A^ — > A^A collision followed by AA^ — > NNcf). No specific model is assumed for 
AA^ — i> NN(f), which is indicated in Fig. Elb) with the serrated line, as we neither did for 
the A^A^ — i> NN(f), but, based on the arguments given above, we will simple assume that 
this process has a cross section about 4.5 times bigger than A^A^ NNcj). 

We evaluate the contribution of the diagram of Fig.|21b) to the nucleon selfenergy in the 
medium. A similar diagrammatic approach was used in [46] to account for corrections to 
ordinary multistep processes but including only nucleons in the intermediate states. How- 
ever, it was rightly concluded there that the equivalent mechanisms with A intermediate 
states should be more important than the mechanisms evaluated there. 

We calculate diagram El b) in two steps. First we evaluate the selfenergy of the incident 
proton due to the generic diagram of Fig. El a) : 

-iE{k) = 3 / -— - — — a-qa- [-q) 



(27r)4 V ; ^ g02_^2_^2 _n^ 
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Figure 3: Diagram for the A^A^ NA reaction for A excitation on the target a). Previous 
mechanism followed by the AA^ NN(f) process b). The long-dashed line represents the 
cut for the calculation of the imaginary part of the selfenergy. 



X 



i{l — n{k — q)) 



+ 



in{k — q) 



ko -qO -EN{k-q)+ie - q^ - E^ik - q) - ie 
where the factor 3 is an isospin factor and IIa represents the pion selfenergy for A/i 



excitation. One can now perform a Wick rotation as it is done in [53], see Fig. 0] The 
integral along the imaginary axis only contributes to the real part of E, and using the 
contour in the complex plane of the variable g° shown in the figure, one picks up only 
the nucleon pole of the 1 — n{k — q) term in Eq. ()18|) . hence g° becomes k^ — E^ik — q). 
One then obtains the imaginary part of the selfenergy, corresponding to the cut shown in 
Fig. El b), as 



TmE{k) = 3 



NN 



d^q 
(2^ 



n(k — q))q' 



2 ' 



(19) 

where XmllA accounts for all possible decay channels of the A (like Nn or Nph or Nph(f)). 

The mechanism of production will pick up from XmllA in the numerator of Eq. ()19|) 
only the process with production, but in the denominator XmU^ will appear accounting 
for all the A decay channels, out of which only the most relevant, Nn, will be kept. 

IIa can be calculated in the same way and a final expression can be found in the 
appendix of [54]. For our purpose it is sufficient to take 

(20) 



nA(g) = g 



4 fU 



NA 



^2 

q p 



il-MA-^SA^A(g)' 



where p is the nuclear density and Sa = (q^ + M)"^ — q^- In Eq. (j20|) and in what follows 
we have neglected the three-momentum, p, of the nucleons in the Fermi sea in the energy 
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k"-EN(k-q) + ie 




V////////////////////(9y//////////////////A 

k''-EN(k-q)-ie 







Re(qO) 



Figure 4: Analytic structure and contour used in the (f integral of Eq. (jl8|) (Wick rota- 
tion). The shadowed regions represent the position of the cuts from the renormalized pion 
propagator. 

denominator since these momenta are small compared to typical values of q and k. Note 
that we are keeping the Mi^jE^Xq) relativistic factors. /^From Eq. we find 

x^n.(,) ^ 1 r^Vrp r#^V — — (21) 



SA - Ma - 



Ma 



For IIa in the denominator of Eq. (fTTH) it is sufficient to use Eq. putting Sa = 
— zr(sA)/2, where r(sA) is the A width for ttA^ decay. However, in the numerator of 
Eq. p9|) . XmllA is given by Eq. (PT|) where Xttt-Sa in the numerator should only account 
for production. In this way, when XmllA is placed in the numerator of Eq. (fT^ it leads 
to the relevant term of production shown in Fig. El b). Thus, we need to evaluate the A 
selfenergy for the process shown in Fig. 

This diagram can be easily evaluated and gives (neglecting in the A four-momentum, 
J9 + g, the momentum p of the nucleon in the Fermi sea): 



-^SA(p + g) = /(l^/ ^HT)(-.T)zW-p')x 
M i{l-n{q-q')) i 



EN{q- q') p° + q^- - E^iq- q') + te p'O^ -p'^ -^l + ie 



(22) 



where U {q' — p') is the Lindhard function [55] with the normalization of the appendix 
of [54], and we have only kept the particle part of the nucleon propagator (the only one 
which contributes to the imaginary part, as we saw when doing the Wick rotation). In 
Eq. (j22j) T is the AA^ —>■ NN(f) amplitude, which we simple assume to be 2.13 times bigger 
than the NN NNcj) one. 
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p+q 



Figure 5: A selfenergy diagram accounting for the AA^ — > NNcj). The long-dashed hne 
represents the cut for the calculation of the imaginary part of the selfenergy. 



Then XmSA can now be obtained again using a Wick rotation, or easier, applying 
Cutkosky rules [56] adjusted for the present problem [57] for the cut shown in Fig. El 

Ea — >■ 2iXmSA , 

U{q' - p') ^ 2ie(g'° - p'^)ImU{q' - p') , 

G{p + q- q') 2iQ{q^ + p° - q'^)ImG{p + q - q') , 

D{p') 2ie{p"^)ImD{p') , 

where G and D are correspondingly nucleon and propagators. These rules lead to the 
expression 



d^q' f d^p' 



(27r)3 J (27r)3 



[l-n{q-q')][l~n{q' -p')]\T\^X 



X 



-(-vr) + 2M-Er,ik-q)-Er,{q- q') - E^iq' - p') - uj^{p')) , (23) 



2uj^{p 

which shows explicitly the (5-function of energy conservation for the process NNN 
NNN(j). In Eq. ()23|) we have used the following approximation: 



ImU{q' - p') ^ -7rp[l - n(g - q')]5{q'^ - p'^ + M - Ei^{q' - p')) (24) 
Altogether we obtain the following expression for IirLL{k) in Eq. (fTU|) : 



\p' 



M 



47r 3 y 



NN 



d^q f d^q' 



(27r)3 J (277)3 



dp'x 



^ Mq^—^—-[l-n{k-q)][l-n{q'-p%l-n{q-q')] 



Ig"'! E^ik - q) EN{q- q') 2a;^(P') 

1 



Ma 



X 



X 



(A;0 _ Ei,{k - q)Y -(P-ml- Yi^{q) 
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^-^^a-^Sa 



F(g)^x 



X 9(1 - cosi9p^-,y) e(A;° + 2M - Er,{k - q) ~ E^{q- q') - u^{p') - M) 



(25) 



with 

cos(^p,-,) 



M'' +p''' + q''' -[k" + 2M - E^ik - q) - E^iq- q') - uj^{v')? 



In Eq. we have aheady exphcitly introduced a form factor F{q) for any ttNN and 
ttNA vertices 

- ml 

^'«) = 

with A = 1 GeV, as usually done. 

Next we want to interpret XmT,{k) in terms of a cross section for production. We 
follow [57, 58] and write: 

r = -2-^Xmi:{k) (26) 

for the probability of production per unit time, keeping the important relativistic factors. 
Then 



dl \k\ 



(27) 



and the cross section for the process, taking into account the initial proton distortion and 
the final distortion, is given by 



(7=1 d^r exp 



Lab 



-2M 



IrnL^pi 



ab 



\Pi 



ab 



f28) 



where TmT,{piab) is calculated with the same expression as in Eq. (j25j) including in addition 
the final (p distortion factor of Eq. (jUj), 



exp 



dl^j^lmll{p^, p{r')) 



(29) 



We simplify the calculations by assuming the to go in the forward direction, p^ = (0, 0, p') , 
which is a good approximation in the nucleus rest frame, where the calculations are done, 
as we have checked. 

Next we briefly describe how to take into account A excitation in the projectile taking 
advantage of the formalism described above. The relevant diagram that we should consider 
is given in Fig. IHl 

The notation for the momenta in Fig. IHl have been chosen to keep maximum analogy with 
the previous mechanism. Hence the same formulae as before can be used adding to 2mT,{k) 
of Eq. for A excitation in the target another term XmE'(fc) to account for A excitation 
on the projectile, which has the same expression as Eq. with the following changes: 
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Figure 6: Diagram for the two-step production with intermediate A excited on the 
projectile. 



^ (A- - g)^, 

F{q)^F{k-q~p), 
1 

- Ej^ik - -q^-ml- UAiq) 

^ ^ ^ . (30) 

{EN{k -q) - My - (k-q)^ -ml-UA{k-q- p) 

Let us note that the last mechanism is similar to the two-step mechanism with only 
nucleons studied above in subsection 12.21 In the latter case, the projectile nucleon loses 
some energy and later on collides with other nucleons to produce the (p. In the mechanism 
of Fig. ini the projectile gets excited to a A and this A collides with other nucleons to 
produce the (p. Since we estimated the cross section for AA^ —>■ NNcf) to be about four 
times bigger than for A^A^ NN(f), we should expect that the mechanism of Fig. should 
be more relevant than the two-step mechanism involving only nucleons. This is indeed the 
case, as we shall see in the results, showing also that the mechanism of A excitation in the 
projectile is more important than that of A excitation in the target in the present reaction. 

It is worth mentioning that this diagrammatic method cannot be directly used for 
the two-step nucleon mechanism, because if a nucleon replaces the A in the diagram of 
Fig. ini this nucleon could be on-shell and would live forever, thus producing formally an 
infinite amount of (j) in its collision with infinite nuclear matter, where the calculations are 
done. This is not the case for the A due to the short lifetime. Hence, in this case the 
infinite matter approach, together with the local density approximation (also used in the 
calculations), can be reliably used. For the nucleon case the explicit consideration of the 
finite size of the nucleus is essential. This is the reason for the different treatment of these 
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two processes in spite of their similarity. These arguments are further elaborated in [59]. 

3 Results and discussion 




1 1 I I I I I I I I I I I I I I I I I I I I I I I I 

20 40 60 80 100 120 140 160 180 200 220 240 

A 



Figure 7: The results of our calculations for proton kinetic energy Tp = 2.7 GeV with only 
one-step process. Different lines correspond to separate contributions from the PS, ISI, 
FSI and total cross section. 

First we want to present results from the model of subsection 12.11 accounting for the 
one-step process. In Figs. [Z|and|Hlwe show the results for R = a a/ {Aa free) for the pro- 
jectile energies Tp = 2.7 and 2.83 GeV respectively. Different lines correspond to separate 
contributions from the phase space (PS), Eq. (0); phase space and initial state interaction, 
Eq. (P) including the distortion factor of Eq. (ISI); phase space and final state inter- 
action, Eq. ((1} including the distortion factor of Eq. (jH)), (FSI); and complete calculation, 
(total), i.e. the simultaneous contribution of all the effects, Eq. 0. We performed calcu- 
lations for the following nuclei: I'^C, 1^0, j^Mg, f^Al, f^Si, f^P, IjS, i°Ca, f^Fe, ^Cu, 

89 V 110/^^ 152 c™ 208 pz, 238tt 
39^ ) 48 '-^"5 62 '-'"h 82 '^J 92 • 
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Figure 8: The same as Fig. [3 but for proton kinetic energy energy Tp = 2.83 GeV. 



By looking at the PS curves of Figs. U\ and |H1 we see that R is larger than unity in 
both cases and much larger for the lower Tp. This is due to the Fermi motion. This 
increase of the nuclear cross section over Aafree is a well known fact in subthreshold 
production in particle- nucleus as well as in nucleus- nucleus collisions [60,61]. The reaction 
threshold is defined for the scattering on a nucleon at rest. Yet, at this threshold energy, 
the Fermi motion of the nucleons makes the reaction possible and the ratio R would 
grow up to infinity as we approach the threshold energy. At subthreshold energies, the 
consideration of large momentum components of the nucleons in the nucleus by means of 
the spectral function, accounting simultaneously for important correlations between energy 
and momentum [38,62], helps to increase the cross section [63], as well as the consideration 
of multi scattering processes [64,65]. All these mechanisms become progressively less and 
less important as we go up to energies above threshold, and as we discuss below, the possible 
uncertainty of our results from neglecting these sources will not modify our conclusions, 
since we are interested in the form of the ^-dependence of the cross section rather than in 
its absolute value. Indeed, the PS curve is practically constant with a good accuracy for 
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all A. The effects discussed above, from the spectral function and multinucleon scattering 
mechanisms, are volume effects, not affected by the distortion of the initial proton or the 
absorption of the final 0. Thus the constancy in A of the corresponding PS calculations 
including these new effects would also hold with a somewhat increased cross section. Now 
if we look at ISI and FSI curves, we see that in both cases there is a sizeable decrease 
of R, particularly for ISI, which shows a stronger A dependence. Although the ISI and 
"total" curves are almost parallel, the absolute values decrease with A and therefore the 
contribution of the FSI becomes more and more important. This significant A dependence 
can be seen in the ratio of these two curves which is shown in Fig. IHl We see that the 




Figure 9: Ratio of the total cross section to ISI with only one-step process. Calculations 
are done for Tp = 2.83 GeV. 

ratio decreases from values around 0.8 for light nuclei to values close to 0.5 for heavy 
nuclei. From this figure we can conclude that in the A-dependence there is indeed valuable 
information concerning the absorption and hence, the width in the medium, which is 
the main conclusion of the present work. 

At this point we want to discuss the effect of taking into account the momentum 
dependence of the nucleon selfenergy. By doing the changes explained in Eqs. ^T^i and 
we find an increase of about 15% in the cross sections (which we would consider as 
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an upper bound), but this percentage increase is remarkably equal for all nuclei such that 
the A dependence of the curves in Figs. [TllHlis preserved. This means that if we normalize 
the cross section to the one of the given nucleus, as we shall do later, the effect discussed 
above does not show up. 

Now we pass to consider the contribution from the two-step processes with nucleon 
and A intermediate states. In Figs. ^1 and E] we plot the ratios of the different two- 
step mechanisms to the one-step process for two different energies as a function of the 
mass number. We see that the most important two-step contribution comes from the 
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Figure 10: Ratio of the nuclear cross section of the different two-step processes over the 
cross section of the one-step process. The different lines represent: dotted line, only nu- 
cleons as intermediate state with AE = 300 MeV (see subsection 12. 2j) : dashed-dotted line, 
only nucleon with AE = 400 MeV; solid line, intermediate A excited on the target (see 
Fig. El; dashed line, intermediate A excited on the projectile (see Fig. E}. The calculations 
are done for a kinetic energy of the beam proton of Tp = 2.7 GeV. 

A excitation in the projectile which is comparable to the one-step mechanism and about 
5 — 10 times bigger than the two-step mechanism involving only nucleons or A excitation 
on the target. Even then, we are concerned about the A dependence, no so much on the 
absolute values of the cross sections, since the absorption effect is reflected in this A 
dependence. 

The relevance of the two step processes in (p production was also put of manifest in 
Ref. [66]. There it was shown that the two step mechanism, in which a pion is produced 
in the intermediate states was dominant below threshold, but small compared to the one 
step process at the energies studied here. We estimate our two step cross section involving 
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Figure 11: Same as Fig. CHlbut for Tp = 2.83 GeV. 



A in the intermediate states to be of the same order as the one step process. We should 
note, however, that the two step processes involving vr or A in the intermediate are not 
the same, in the sense that the serrated line in Figs. 01 and IHl does not necessary stand 
for an on shell pion. Actually, as shown in [46], having in this line intermediate p states 
leads to a much larger two step contribution. We accept uncertanties from the two step 
mechanism, but then we look for one observable that minimizes these uncertanties and 
this is given by the A dependence of the cross section. Indeed, the fact that the curves in 
Figs. ^1 and ^2 are almost flat as a function of A indicates that the A dependence of the 
sum of all mechanisms has essentially the same A dependence as the one-step mechanism 
alone. To see this more clearly, we show in Fig. the ratio R{^X) normalized to R{}'^C) 
for the one-step and one- plus two-step mechanisms. What we see in the figure is that 
this normalized R changes very little when including the two-step mechanisms for both 
the Tp considered. Note that for the energy closer to the threshold (Tp = 2.7 GeV) the 
changes due to the two-step contributions are smaller. Therefore we conclude that the A 
dependence obtained in the present work is reliable and the calculations clearly show that 
proton induced (p production in nuclei at energies just above threshold can indeed be used 
to get information on the (f) width in the medium. 

In order to see which is the experimental precision needed to get a definite information 
on the (f) width in the medium, we have performed the same calculations assuming widths 
in the medium to be one half or twice the width used so far [8, 14]. In the calculations 
this is implemented by multiplying the argument of the exponential in Eq. (jU)) by 1/2 
or 2 respectively. In Fig. [THlwe show the results of these calculations for Tp = 2.83 MeV 
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Figure 12: Ratio of the nuclear cross section normalized to for two different incident 
proton kinetic energies, Tp, including or not the two-step mechanisms. The two-step process 
with nucleon intermediate states has been evaluated with AE = 400 MeV. 

(without the inclusion of the two-step processes). Comparing Figs. ^1 and ^1 we clearly see 
that the uncertainties due to the second step mechanism are far smaller than the differences 
in the results obtained by using these different widths. The three curves shown there for 
the width of [8, 14], half this width, and double, should serve to get a fair answer about 
the (f) width in the medium by comparing with experimental results. Comparing Figs. IT^ 
and one can see that the uncertainties one might have from the approximate knowledge 
of the two-step processes still would allow us to be sensitive to the value of the width in 
the medium to the level of 25% of the (f) width we have used. 

Now we address another question, which has to do with (p production in N ^ Z 
nuclei. The calculations of the (p production nuclear cross sections are done in symmetric 
nuclear matter. Hence in order to calculate the relative (p production cross section R = 
a a/ {Aa free), wc implicitly took a total free elementary production cross section a free = 
{o-pn,(t> + o'pp,(t>)/'2, therefore in a strict sense our model is valid only for the nuclei with equal 
amount of p and n, i.e. up to i^Ca in our calculation. Since the averaged |Tp used for 
elementary production cancels in the numerator and denominator of R, our model can 
also be considered for any other series of nuclei with the same ratio of Z/N. Experimentally 
we have poor knowledge about these elementary cross sections: there is experimental 
data [35] for the pp ^ pp (p reaction at one energy, close to the one used in our work, and 
nothing for pn ^ pn (p. Some models (see for example [45]) tell us that <7pn,(f)/o'pp^(f, ~ 5 
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Figure 13: Ratio of the nuclear cross section normalized to for Tp = 2.83 GeV multi- 
plying the width in the medium by different factors. 

for our energies. Nevertheless, our results can still be used to compare with experiment 
if one takes for a free the isospin weighted combination {Napn,<f, + Zapp^^)/A. Should we 
know these elementary cross sections we could compare with the experimental nuclear cross 
sections ioi N ^ Z nuclei multiplying the presents results for R by [Napn^^ + Zapp^^)/A. 
Note, however, that, even with Opn^^j Opp^^ = 5, the ratio 

{Napn,^ + Zapp^^)/A ^^^^ 

for a very asymmetric nucleus like ^ssjj jg j^g^ 1.15, a small correction compared to the 
effects from (p absorption in this nucleus. However, this correction is not so small if we 
consider that the difference in Fig. for A ~ 240 for the full width or twice this width 
is only of the order of 33% and between the full width and half this width is of the order 
of 25%. Hence, in order to determine the medium width with a precision of better than 
50% from heavy nuclei the use of this isospin correction is important. 
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On the other hand, if we take a nucleus hke Ca, which is isospin symmetric, the ratios 
of the curves with 2r, F and O.SF to the one with OF are 1.50, 1.35 and 1.20 respectively. 
This gives us an idea of the precision one can get for F^ given a certain precision in the 
experimental results. In our opinion, a way to get a high precision on this experimental 
ratio would be to make a fit to a data set for several approximately symmetric nuclei and 
obtain the ratios from the fitted curve. 



4 Conclusions 

We have performed calculations of relative cross sections for (f) production in nuclei in 
proton nucleus collisions with the aim to obtain information on the (p width in the nuclear 
medium. For this purpose we explored the A dependence of the cross section which is 
tied to the absorption of the (p iii its way out of the nucleus. In the absence of initial 
proton and (f) distortions the cross sections obtained are practically proportional to A. 
Sizeable diversions from this linear dependence come when both distortions are considered. 
Although the initial state interaction was more effective reducing the cross section, even 
then we found sizeable changes due to absorption which result in an extra reduction of 
the cross section in about a factor of two in large nuclei. These predicted changes are large 
enough such that devoted experiments can obtain relevant information on the (p width in 
the medium. Since the A dependence of the cross sections, and not so much the absolute 
values, is important to learn about (p absorption, we present cross sections for heavy nuclei 
normalized to a light one, which should be specially suited for comparison with future 
experiments. We have also seen that in order to extract the optimum information on the 
(p width it would be useful to have data on (p production on neutron targets, for what 
experiments on the deuteron would also be most welcome. The calculations have been 
done at an energy just above threshold, which is accessible in the COSY facility at Juelich, 
and for which our theoretical treatment of the initial state interaction is easy and reliable. 
The results obtained in the present work clearly show that the modification of the (p width 
in the nuclear medium has sizeable effects in this reaction to the point that the actual 
experimental implementation of the reaction should provide a measure of the strength of 
the medium (p width. 
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